Br | 


干旱 区 研究 


ARID ZONE RESEARCH 


doi:10. 13866/j. azr.2019.02.23 


我 国 荒漠 植被 生产 力 动 态 及 其 与 水 热 因子 的 关系 、 


NSE, BB, R 


A, 杜 自强 *， 张 ” 红 : 


(1. 山西 大 学 环境 与 资源 学 院 , 山 西 太原 030006; 2. 山西 大 学 黄土 高 原 研究 所 ,山西 太原 030006; 


3. 西北 农林 科技 大 学 资源 与 环境 学 院 , 陕 西 ve 


712100) 


摘 要 : 为 了 探讨 近 30 a 来 我 国 干旱 区 荒漠 植被 的 净 初 级 生产 力 (NPP) 及 其 与 水 热 因 子 相关 性 随时 间 的 变化 , 运 
用 CASA( Carnegie-Ames-Stanford approach ) 模型 估算 我 国 荒漠 植被 1982—2015 年 生长 季 的 NPP, 并 运用 线性 回归 和 
GIS 空间 分 析 方 法 分 析 了 NPP 的 时 空 变化 特征 ,利用 滑动 相关 系数 分 析 了 荒漠 植被 NPP 与 水 热 因 子 的 关系 。 结 


RRI: 单位 面积 NPP 均值 为 42 gm? 


“a NPP 整体 水 平 较 低 。 空 间 上 时 西北 部 、 东 部 边缘 较 高 ,中 部 、 南 
部 和 中 东部 较 低 的 分 布 特征 。@) 荒漠 植被 NPP 年 均 总 量 为 5.783 x 10" g 


a-!。 从 落 漠 植被 NPP 的 年 际 变化 来 


看 ,1982 一 2015 年 中 国 荒漠 植被 NPP 总 量 以 1.64 x10”g. (10a) 的 线性 速率 ( 忆 = 0.054) 上 升 ,荒漠 植被 生长 状 
况 总 体 上 不 断 改善 ,但 总 量 趋势 呈现 阶段 性 变化 ,1982 一 1993 年 荒漠 植 被 NPP 总 量 呈 极 显 著 增 长 态势 (1.25 x 10” 


ga ,PP<0.01);1993 一 2006 年 NPP 总 量 呈 极 显 著 降 低 态 势 ( -6.42 x 10" g+ am! P <0. 01) ;2006—2015 年 


NPP 总 量 缓慢 增长 (1.70 x10" g .a”,P>0.05)。 从 空间 变化 来 看 ,47.65% 的 荒漠 植被 NPP 呈 增 加 态势 ,主要 分 
布 在 阿拉 善 高 原 、 天 山北 区 塔里木 盆地 西部 边缘 、 柴 达 木 盆地 的 东南 边缘 阿尔 金 山南 甘 和 昆仑 山脉 。@) 从 荒漠 
植被 NPP 与 各 气候 因子 之 间 的 相关 关系 随时 间 的 变化 来 看 ,NPP 与 气温 的 滑动 相关 系数 随时 间 的 变化 保持 为 负 
相关 ,与 降水 .干燥 度 的 滑动 相关 系数 保持 为 正 相 关 ,与 太阳 总 辐射 的 滑动 相关 系数 随时 间 变 化 并 未 表现 出 显著 的 
变化 趋势 。 总 体 上 ,荒漠 植被 与 水 热 因 子 的 相关 关系 在 研究 时 段 均 有 进一步 减弱 的 态势 , 即 荒漠 植被 NPP 对 气候 


因子 的 变化 愈 来 愈 不 敏感 。 


关键 词 : 净 初级 生产 力 ; MRE; 气候 因子 ; 滑动 相关 ; 水 热 因 子 


净 初 级 生产 力 (net primary productivity ,NPP) 是 
指 绿色 植物 单位 时 间 .单位 面积 所 累积 的 有 机 干 物 
质 总 量 ,是 由 光合 作用 所 产生 的 有 机 质 总 量 中 扣除 
自 养 呼吸 后 的 剩余 部 分 '" ,不仅 能 反映 自然 植被 的 
生产 能 力 和 陆地 生态 系统 的 质量 ,而 且 在 测量 陆地 
生态 系统 可 持续 发 展 能 力 方面 也 起 着 重要 作用 "”。 

荒漠 植被 是 陆地 生态 系统 的 重要 组 成 部 分 , E 
要 分 布 在 我 国 的 西北 干旱 地 区 ,是 在 极端 干旱 和 恶 
劣 的 生境 下 发 育 的 植被 类 型 ,具有 发 达 的 根系 和 独 
特 的 水 分 生理 特征 ,抗旱 能 力 较 强 ,也 能 较 快 的 响应 
外 界 环境 的 变化  。 前 人 有 关 荒 漠 植 被 生长 状况 
与 气候 因子 关系 的 研究 表明 ,我 国 西北 干旱 地 区 过 
去 30 年 气候 向 暧 湿 化 转型 ,植被 主要 旦 改善 趋势 ， 
降水 量 和 莹 散发 对 植被 生长 影响 较 大 ” ; 吴 晓 全 
等 ”模拟 的 天 山 荒漠 NPP 为 55.78 gm a， 
NPP 总 量 波 动 趋 于 平缓 , 且 若 漠 植 被 NPP 主要 受降 
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雨量 的 影响 ; 叶 静 芸 等 '" 采 用 RVI 线性 模型 估算 我 
国 乌 兰 布 和 沙漠 东北 缘 的 鞠 漠 -绿洲 过 渡 带 内 芳 漠 
植被 的 单位 面积 地 上 生物 量 为 90.73 g mm 一; 李 芳 
等 ”研究 表明 ,黑河 中 游 荒漠 植被 生产 力 对 降水 响 
应 具有 兆 后 期 特征 ;Cui 等 ”认为 ,阿拉 善 和 河西 走 
廊 植 被 生长 状况 表现 出 退化 趋势 ,内 蒙古 高 原 中 部 ， 
西部 大 部 分 地 区 植被 则 呈现 增长 趋势 ; 张 凯 等 ” 认 
为 ,阿拉 善 地 区 东部 车 漠 植 被 有 所 改善 ,但 中 西部 地 
区 出 现 退 化 现象 ,降水 量 对 荒漠 植被 的 控制 作用 较 
强 ; 也 有 研究 认为 ,西北 干旱 区 植被 生长 状况 与 降 
水 ,温度 的 相关 性 都 不 高 "… ” 。 以 上 的 研究 大 多 集中 
于 西北 干旱 区 菏 一 区 域 的 总 漠 植 被 的 生长 状况 , 且 
进行 植被 与 气候 因子 关系 分 析 时 关注 的 是 整体 的 相 
关 性 ,不 能 反映 植被 与 气候 因子 相关 关系 的 时 间 变 
化 。 因 此 ,在 西北 地 区 暧 湿 化 背景 下 ,进行 长 时 间 序 
列 ,区 域 尺 度 上 芝 漠 植被 的 NPP 格局 及 其 与 气候 因 
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子 相 关 性 随时 间 变 化 的 研究 十 分 必要 ,这 将 有 助 于 
全 面 认识 植 被 一 气候 相互 关系 ,揭示 气候 变 暖 给 陆 
地 生态 系统 带 来 的 影响 ” ,为 荒漠 生态 环境 调节 与 
生态 恢复 提供 科学 依据 。 


1 材料 与 研究 方法 


1.1 研究 区 概况 

荒漠 植被 分 布 区 (图 1) 是 通过 ArcGIS 软件 从 
我 国 1:4 000 000《 中 国 植被 图 》 中 提取 出 来 的 ,主要 
分 布 于 我 国 的 西北 干旱 区 , 介 于 73° ~ 115° E 和 
35° ~50°N, 为 年 平均 降水 量 小 于 200 mm 且 地 理 位 
置 在 贺兰山 以 西 ,总 面积 约 2.5 x10 km 的 广大 区 
域 。 区 域内 山 盆 相 间 ,西高 东 低 , 南 高 北 低 。 和 常年 受 
大 陆 性 气候 的 影响 ,降雨 稀少 , 蕉 发 强烈 ,是 世界 上 
最 为 干旱 的 地 区 之 一 。 植 被 主要 是 以 早生 、 超 早生 
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精度 。 本 文 进一步 采用 最 大 值 合成 法 消除 了 一 部 分 
云 大 气 的 干扰 ,并 减少 月 内 物候 循环 的 影响 "” ,得 
到 每 月 的 NDVI 数据 集 。 

荒漠 植被 包括 温带 矮 半 汐 木 启 漠 Ui VT ER 
EPR TEAR IR Ty HEE ANE HE AR ie UTS AE FRAN 
Fre Ve A Hla HF Fe FE FA IEEE ARTE GREET 
中 国 科 学 院 植物 研究 所 编制 的 1:4 000 000《 中 国 植 
被 图 》, 通 过 ArcGIS 软件 提取 。 

本 文 所 需 的 大 气 驱 动 数据 包括 太阳 总 辐射 月 
平均 气温 、 月 降水 量 .干燥 度 。 数 据 来 源 于 中 国 气象 
科学 数据 共享 服务 网 提供 的 中 国 地 面 气候 资料 月 值 
数据 集 。 干 燥 度 采用 基于 气温 与 降水 数据 计算 的 适 
于 较 大 尺度 研究 广泛 采用 的 de Martonne 干燥 度 指 


数 "”。 在 此 基础 上 ,通过 反 距离 权重 法 "对 各 个 
气象 因子 进行 空间 插值 ,生成 与 NDVI 数据 具有 相 


的 灌木 . 半 灌 木 及 小 乔木 为 建 群 种 构成 的 稀疏 植被 ， 


群落 中 植物 种 类 稀少 ,区 系 成 分 单调 ;该 区 域 又 因 地 
处 中 亚 、 西 但 利 亚 .蒙古 、 ee nn ed aE 
自然 地 理 条 件 在 历史 上 又 几经 变迁 ,为 各 类 植物 区 
系 成 分 的 接触 .混合 和 特 化 提供 THAR RE, EHK 
系 特征 具有 明显 的 特殊 性 5 。 
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图 1 荒漠 植被 分 布 区 


Fig.1 Distribution areas of desert vegetation 


1.2 数据 与 方法 

1.2.1 数据 来 源 及 处 理 ”遥感 数据 是 美国 国家 航 

天 航空 局 全 球 监测 与 模型 研究 组 (CIMMS ) 提供 的 

NDVI3g 数据 。 时 间 序 列 为 1982 一 2015 年 ,空间 分 

辩 率 约 为 8 km x8 km ,时 间 分 辩 率 为 15 d。 该 数据 

集 消 除了 太阳 高 度 角 传感器 灵敏 度 随 时 间 变 化 等 
的 影响 ,并 结合 交叉 辐射 定 标的 方法 ,增强 了 数据 的 


同 空间 分 辨 率 和 投影 方式 的 气象 栅 格 数据 。 
1.2.2 研究 方法 

(1) CASA 模型 

Æ x F] H CASA ( Carnegie-Ames-Stanford 
approach ) 模型 来 估算 我 国 荒漠 植被 1982—2015 年 
的 生长 季 NPP? 。 该 模型 主要 涉及 2 个 变量 ,分 别 
为 植被 吸收 的 光合 有 效 辐 射 (4PAR) 和光 能 利用 率 
(e)s AMF: 

NPP(x,t) =APAR(x,t) xe(x,t) (1) 
INP APAR (x,t) RIR tA MRC x 处 吸收 的 光合 有 
效 辐射 (MJ m?) ze (x,t) RIR t HORIE x 处 实 
际 光 能 利用 率 (g . MJ) 。 

光合 有 效 辐射 计算 公式 如 下 : 

APAR(x ,t) =SOL(x,t) x FPAR(x,t) x0.5 (2) 
rf SOL (ast) EAR t 月 份 像 元 x 处 的 太阳 总 辆 身 
量 (MI + m?) ;FPAR(%,t) 表示 植被 层 对 入射 光合 
有 效 辐射 的 吸收 比例 ;0.5 为 植被 能 利用 的 太阳 有 
效 辐射 (波长 0.4 ~0.7 pm) 占 太阳 总 辐射 的 比例 。 

FPAR 由 NDVI 和 植被 类 型 来 表示 ,并 不 超过 
0.95, 


SR( x,t) -SR pin 


oc: a SR max = SR min : 


0.95 (3) 


WP: SR (x,t) RIR t HRI x 处 的 比值 指数 ; 
SR pin BUEN 1.08, SR na HI KD FEL BRE FY ASK, BL 
值 范围 在 4.14 ~6.17, SR(x,t) H NDVI (x,t) 3R 


4 


得 : 
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1 + NDVI(x ,1) 
1 — NDVI(x,t) 
光 能 利用 率 (e) 受 温度 和 水 分 条 件 的 影响 , 它 指 的 
是 植被 将 吸收 到 的 光合 有 效 辐射 (PAR) 转 化 为 有 机 
碳 的 效率 ,计算 公式 如 下 : 

e(x,t) =Ta (x,t) XToCxt) xW, (x,t) X Enx (5) 


RP: Ta (%,1) ,Ta (%,1) 表 示 温度 对 光 能 利用 率 的 
W, 表示 水 分 对 光 能 利用 率 的 影响 ;s, 表 示 理 
想 状态 下 的 最 大 光 能 利用 率 (g + MI“) 。 
温 胁迫 因子 的 估算 :Ts Cest) 和 Ta (x,t) BIR 
温度 对 光 能 利用 率 的 影响 。 
Ta (x,t) =0.8 +0. 02 xT, (x) -0.0005x 
[Ti (x) ]? (6) 


式 中 :7,(x) 表 示 某 一 区 域 一 年 内 NDVI 值 达 到 最 
高 时 的 当月 平均 气温 (%C ) , 当 某 月 的 平均 气温 小 于 
或 等 于 -10 CC 时 ,7T,,(x) 取 0。 


SR(x,t) = (4) 


1. 184 


椿 2 认为 ,以 10 ~20 a 的 滑动 年 数 效果 最 好 。 因 
此 ,为 了 确定 荒漠 植被 NPP 与 气候 之 间 关 系 的 时 间 
变化 ,笔者 选择 17 a 滑动 步 长 ,分 别 计算 了 生长 季 
NPP 和 气温 、 降 水 、 太 阳 总 辐射 .干燥 度 之 间 的 相关 
系数 ,得 到 1982—1998 年 、1983 一 1999 年 ，……. , 
1999—2015 年 各 时 段 的 相关 系数 序列 , 据 此 来 分 析 
我 国 荒漠 植被 NPP 与 气候 因子 相关 系数 随时 间 变 
化 的 稳定 性 特征 。P <0. 05 表示 相关 性 显著 ;P < 
0.01 表 示 相 关 性 极 显著 。 

(3) 线性 回归 分 析 

采用 一 元 线性 回归 拟 合 荒漠 植被 NPP SAR 
子 和 滑动 相关 系数 的 变化 速率 ,通过 i 检验 来 检验 
变化 速率 (线性 回归 系数 ) 的 显著 性 ")。 


2 结果 与 分 析 


2.1 匾 漠 植被 NPP 的 时 空 变 化 
1982—2015 年 中 国 荒漠 植被 单位 面积 年 均 
NPP 为 42 g m”? , 呈 西北 部 .东部 边缘 较 高 ,中 部 、 


Talt) =T ep 0.2 x Tp (2) -10 -TtT 


1 
{1 +exp[0.3 x ( -T,,,(%) -10 -T(x,t))]} (T) 


如 果 某 月 的 平均 气温 7(x,t) 高 于 最 适 气温 
Talx) =10 或 低 于 T(x) =13 CRF, WIA AY 
7,,(x,t) 值 是 月 平均 气温 7 了 (x,t) H Toa (x) 时 
T(x,t) 的 1/2。 

水 分 胁迫 因子 的 估算 : 取 (x,i) 反 映 的 是 水 分 
对 植物 光 能 利用 率 的 影响 , 随 着 环境 中 有 效 水 分 的 
增加 ,到 (xz , 划 逐 渐 增 大 , 它 的 取 值 范围 为 0.5( 在 极 
端 干旱 条 件 下 ) 到 1 (非常 湿润 条 件 下 ) 。 

W,(«,t) =0.5 40.5 xE(x,t)/E,(%,t) (8) 


RP E(x ,t) ers KS AE (mm) ,可 根据 周 
T EE AK, Și EW, AY K aS Di EI ORO”? ; 
E,, (x,t) RIR KEES BCE, ATH YE Bouchet 提 
出 的 互补 关系 求 取 。 

最 大 光 能 利用 率 的 确定 :根据 误差 最 小 原则 , 利 
用 中 国 的 NPP 实测 数据 ,模拟 了 各 植被 类 型 的 最 大 
光 能 利用 率 ,该 研究 成 果 在 估算 中 国 植被 NPP 时 得 
到 广泛 应 用 ” o BX su 的 取 值 也 参照 这 一 成 果 ， 
其 中 荒漠 8na N 0.542 g > MJ。 

(2) 滑动 相关 系数 

滑动 相关 系数 用 来 检验 相关 是 否 稳定 。 林 学 


南部 和 中 东部 较 低 的 空间 分 布 特征 (图 2)。 从 单位 
面积 年 均 NPP 的 分 布 状况 来 看 ,NPP 低 于 50 g - 

m? 的 区 域 占 中 国 荒漠 植被 总 面积 的 79. 65% ,分 布 
于 除 天 山 准噶尔 盆地 ABE LL HRS LL BBR Be OT 
高 原 之 外 的 大 部 分 地 区 , 介 于 50 ~ 100 g- mH 
11.05% ,集中 分 布 在 准 格 尔 倪 地 、 贺 兰 山 、 鄂 尔 多 斯 
高 原 西北 部 。NPP 介 于 0 ~ 100 gm 一 的 荒漠 植被 
面积 占 到 了 我 国 荒漠 植被 总 面积 的 90. 70% ,表明 
我 国 荒漠 植被 生产 力 水 平 较 低 。 

1982 一 2015 年 中 国 荒漠 植被 NPP 年 均 总 量 》 
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图 2 1982—2015 年 荒漠 植被 单位 面积 NPP 分 布 
Fig.2 NPP Distribution of desert vegetation per unit 
area from 1982 to 2015 
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图 3 1982—2015 年 荒漠 植被 NPP 总 量 的 年 际 变化 
Fig.3 Interannual change of total NPP of desert vegetation 
from 1982 to 2015 
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少 (P<0. 05) 和 不 显著 减少 的 面积 分 别 占 荒漠 植被 
总 面积 的 18.21% .9. 12% 和 25. 02% ,分 布 在 昆仑 
山 和 阿尔 金山 北 厅 、 塔 里 木 盆地 中 部 .吐鲁番 盆地 、 
准噶尔 盆地 ALC LL 河西 走廊 以 北 地 区 (图 4) 。 
2.2 荒漠 植被 NPP 与 水 热 因 子 相关 关系 的 分 析 

影响 植物 生长 的 气候 因子 有 很 多 , 本文 选择 气 
温 、 降 水 量 、 太 阳 总 辐射 以 及 表征 气候 综合 因子 的 干 
燥 度 指数 来 分 析 荡 漠 植被 NPP 与 气候 变化 的 滑动 
相关 关系 。1982 一 2015 年 荡 漠 植被 覆盖 区 气温 、 太 
阳 总 辐射 分 别 以 0.491 C + (10a)-!'(P<0.01) 和 
60.461 MJ . (10a)-!' (P<0.01) 的 速率 极 显著 增 
加 ,降水 以 5.576 mm + (10a) ~! (P <0.05) 的 速率 
显著 增加 (图 5)。 干燥 度 分 布 在 0.63 ~1.04 之 间 
且 都 小 于 10 ,表明 匾 漠 植被 分 布 区 为 严重 干旱 地 
区 ,并 且 干 燥 度 以 0.02 . (10a) (已 >0. 05) 速 率 
波动 增加 ,这 与 西北 干旱 区 气候 暖 湿 化 结论 相 一 致 。 
本 文 以 17 a 为 滑动 步 长 分 别 计 算 了 1982 一 2015 年 
荒漠 植被 NPP 与 气温 .降水 量 .太阳 总 辐射 .干燥 度 
之 间 的 滑动 相关 系数 ,以 此 来 分 析 荒 漠 植 被 NPP 与 
各 气候 因子 之 间 的 关系 。 

荒漠 植被 NPP 与 气温 的 滑动 相关 系数 (图 6) 
随时 间 出 现 了 相关 符号 的 改变 。1982 一 1998 年 和 


— 1986—2002 年 期 间 ,两 者 呈正 相关 ,但 都 没有 通过 
N 图 4 1982—2015 年 荒漠 植被 NPP 变化 趋势 显著 性 检验 P < 0 05 的 显著 性 检验 1987 2003 年 之 后 两 者 
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> Fig.4 Significance test on the change trend of NPP of desert 


vegetation from 1982 to 2015 


5.78 x10” g-a',W1.64x10" g + (10a) '(P= 
0.054 ) 的 线性 速率 增加 。NPP 年 总 量 呈 现 阶段 性 
变化 ,1993 年 之 前 荒漠 植被 的 NPP 总 量 以 1.25 x 
10”g. a 的 线性 速率 极 显著 增长 (P <0. 01 ) ; 
1993—2006 年 NPP 总 量 却 以 6.42 x10" g- a! AY 
线性 速率 极 显 著 降低 (P <0. 01);2006 年 以 后 的 
NPP 总 量 叉 开始 缓慢 增长 (P >0. 05 ) (图 3)。 逐 像 
元 来 看 ,年 际 荒漠 植被 NPP 以 0.09g.m”*.a! 的 
速率 增加 ,但 变化 速率 表现 出 明显 的 空间 异 质 性 。 
NPP 呈 增 加 态势 的 区 域 占 荒漠 植被 总 面积 的 
47.65% ,其 中 NPP 时 极 显 著 增 加 (P<0.01) .显著 
(P<0.05) 增 加 和 不 显著 增加 的 面积 分 别 占 荡 漠 植 
被 总 面积 的 13.38% .8.17% 和 26. 10% ,主要 分 布 
在 阿拉 善 高 原 、 天 山北 基 塔里木 例 地 西部 边缘 、 柴 
达 木 盆地 的 东南 边缘 、 阿 尔 金 山南 入 和 昆仑 山脉 。 
NPP 呈 降 低 态 势 的 区 域 占 荒漠 植被 总 面积 的 
52.35% ,其 中 ,NPP 极 显著 减少 (P<0. 01) .显著 减 


iat 的 相关 性 发 生 了 转变 ,由 正 相 关 转 变 为 负 相 关 , 且 负 


相关 随时 间 变 化 逐渐 加 强 ;1992 一 2008 年 和 1993 一 
2009 年 的 负 相 关 达 到 了 P =0. 05 的 显著 性 水 平 , 呈 
显著 负 相 关 ;1993 一 2009 年 之 后 负 相 关 关 系 逐 渐 减 
弱 ,并 有 进一步 减弱 的 趋势 。 

荒漠 植被 NPP 与 降水 量 ( 图 6) 在 各 个 年 段 均 
呈正 相关 ,但 相关 系数 的 趋势 显现 出 阶段 性 的 差异 。 
1982—1998 年 至 1987 一 2003 年 相关 系数 持续 降 
低 ,但 没有 出 现 相 关 符 号 的 转变 。1987 一 2003 年 至 
1996 一 2012 年 NPP 与 降水 量 的 正 相 关 关 系 逐 渐 增 
强 , 并 在 1996 一 2012 年 达到 了 显著 正 相 关 。1996 一 
2012 年 之 后 正 相 关 关 系 急剧 减弱 ,有 可 能 会 进一步 
的 减弱 。 

车 漠 植 被 NPP 与 干燥 度 ( 图 6) 在 各 个 年 段 均 
呈正 相关 ,1982 一 1998 年 和 1983—1999 年 正 相 关 
增强 ,1983 一 1999 年 至 1987 一 2003 年 NPP 与 干燥 
度 的 正 相 关 逐 渐 减 弱 。1987 一 2003 年 之 后 相关 系 
数 的 变化 趋势 与 NPP 和 降水 相关 系数 的 变化 趋势 
一 致 。1987 一 2003 年 至 1996—2012 年 的 相关 系数 
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刘 雪 佳 等 :我 国 蕊 漠 植 被 生产 力 动态 及 其 与 水 热 因 子 的 关系 
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图 5 气候 因子 变化 趋势 


Fig.5 Change trends of climatic factors 
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图 6 17a 滑动 相关 系数 


Fig.6 Seventeen-year sliding correlation coefficients 


普遍 比 NPP 与 降水 的 相关 系数 高 , 且 相 关 性 逐渐 增 
强 , 到 1996 一 2012 年 达到 了 显著 正 相 关 。1996 一 
2012 年 之 后 正 相 关 关 系 减弱 ,并 可 能 会 进一步 的 减 
弱 。 

荒漠 植被 NPP 与 太阳 总 辐射 的 滑动 相关 系数 
(图 6) 随 时 间 出 现 了 相关 符号 的 改变 。1982 一 1998 
年 至 1983—1999 年 期 间 ,两 者 正 相 关 关 系 增强 ,都 
通过 了 P <0. 05 的 显著 性 检验 。1983 一 1999 年 至 
1993 一 2009 年 ,两 者 的 正 相 关 关 系 逐 渐 减 弱 , 但 相 
关系 数 在 1984—2000 年 .1987 一 2003 年 和 1989 一 
2005 年 达到 了 P = 0. 05 显著 正 相 关 。1993 一 2009 
年 之 后 ,两 者 的 相关 系数 由 正 相关 转变 为 负 相 关 , 直 
到 1995—2011 年 之 后 相关 系数 转变 为 正 相 关 并 逐 
渐 增 强 ,但 到 1998 一 2014 年 之 后 两 者 的 正 相 关 关 系 
逐渐 减弱 ,并 有 进一步 减弱 的 趋势 。 

总 的 来 说 ,荒漠 植被 NPP 与 气温 的 滑动 相关 系 


数 随时 间 的 变化 主要 保持 为 负 相 关 , 充 漠 植被 NPP 
与 降水 .干燥 度 的 滑动 相关 系数 随时 间 的 变化 保持 
为 正 相关 ,而 荒漠 植被 NPP 与 太阳 辐射 的 相关 系数 
随时 间 未 表现 出 明显 的 相关 性 ,但 四 者 的 相关 关系 
均 有 进一步 减弱 的 趋势 , 即 证 明 荒 漠 植 被 NPP 对 气 
修 因子 的 变化 愈 来 愈 不 敏感 。 


3 ”讨论 与 结论 


(1) 1982—2015 年 ,中 国 荒漠 植被 单位 面积 
NPP 均值 为 42 g+m?+a', AIC 2000—2015 
年 的 估算 结果 与 MODIS 同期 的 500 m 分 辩 率 NPP 
产品 (MOD17A3H ) 进行 了 对 比 , 发 现 MOD17A3H 
在 大 部 分 荒漠 植被 覆盖 较 少 的 地 区 出 现 缺 值 现象 ， 
如 塔克拉玛干 沙漠 .昆仑 山脉 .准噶尔 盆地 BIR TS AN 
盆地 周边 及 阿拉 善 高 原 的 大 部 分 地 区 扣 。 在 与 本 
文 研究 范围 重合 的 3 220 个 像 元 中 ,两 者 的 相关 系 


FOR 区 


BON 0.57 ,达到 了 极 显著 相关 ,说 明 MODIS NPP 产 
品 与 笔者 估算 的 NPP 结果 非常 吻合 。Wang 等 中 人 
算 的 我 国 干旱 半 干 旱地 区 2000—2010 年 的 年 均 
NPP 为 159.31 gm , 昌 在 这 10a 内 NPP 呈 上 升 
趋势 ,但 由 于 研究 范围 包括 水 热 条 件 较 好 的 内 蒙古 
东部 地 区 , 所 以 估算 的 NPP 年 均值 较 高 。 焦 伟 
等 基于 2000 一 2014 年 的 MODIS NDVI 数据 模拟 
的 我 国 西北 干旱 区 荒漠 植被 的 年 均 NPP 为 51.1 
gm ,Pan 等 2 模拟 的 我 国 荒漠 植被 的 年 均 NPP 
为 56.3 gm 一 ,本 文 研究 结果 与 之 相近 ,可 见 ,本 
研究 结果 能 反映 我 国 荒漠 植被 的 NPP 

(2) 过 去 的 34 a 来 我 国 荒漠 植被 NPP 年 平均 
总 量 为 5.782 x 10° g, NPP 总 量 呈 现 出 增长 的 态 
势 ,这 与 北半球 陆地 植被 NPP 在 过 去 几 十 年 有 所 增 
Ja? 的 趋势 相 一 致 ,表明 我 国 荒漠 植被 在 固 碳 方面 
发 挥 着 越 来 越 大 的 作用 。 整 个 时 期 荒漠 植被 NPP 
年 总 量 呈 现 阶 段 性 变化 ,1982 一 1993 年 落 漠 植被 
NPP 总 量 呈 极 显著 增长 态势 (1.25 x 10" ge ant, 
P <0. 01) ;1993 一 2006 年 NPP 总 量 时 极 显著 降低 
态势 ( -6.4 x10" g + a7! P <0.01);2006 一 2015 
年 NPP 总 量 缓慢 增长 (1.7 x10"g. a!',P> 
0. 05 ) 。 通 过 分 析 这 3 个 阶段 的 气温 、 降 水 量 .干燥 
度 和 太阳 总 辐射 的 变化 趋势 ,发 现 NPP 的 变化 趋势 
与 同期 气温 呈 负 相关 ,与 同期 的 降水 ,干燥 度 和 太阳 
辐射 呈正 相关 。 有 研究 表明 ,我国 在 1988 年 .2000 
年 .2003 年 .2010 年 也 都 发 生 过 不 同 程度 的 干旱 , 当 
时 我 国 的 NPP 明显 减少 。 然 而 ,我 们 并 没有 观 
察 到 在 这 几 个 年 份 里 荒漠 植被 NPP 发 生 明 显 下 降 ， 
这 表明 干旱 可 能 并 没有 对 荒漠 植被 的 生产 力 造成 显 
著 的 影响 。 

(3) 我 国 荒漠 植被 NPP 整体 水 平 较 低 , 介 于 
0~100¢-m~?-a 的 面积 占 到 了 90.7% ,这 主要 
是 由 于 我 国 西北 地 区 气候 极为 干燥 ,植被 的 生产 力 
较 低 ” 。 单 位 面积 NPP 大 于 100 g +m? + a! 的 
地 区 主要 分 布 在 天 山 和 祁连山 地 区 ,这 主要 受 局 部 
地 形 和 海拔 的 影响 ,有 较 多 的 降水 二 , 且 降 水 与 
NPP 保持 为 正 相 关 ,降水 的 增加 能 促进 车 漠 植物 的 
生长 。 而 NPP 空间 变化 的 异 质 性 则 与 气温 、 降 水量 
的 地 区 分 布 不 均 和 空间 年 际 变化 速率 不 同 有 关 。 

(4) 笔者 研究 结果 显示 ,荒漠 植被 NPP 与 气 
温 、 降 水 太阳 总 辐射 .干燥 度 之 间 的 相关 系数 具有 
时 间 变 异性 ,这 种 变化 表明 荒漠 植被 活动 对 气候 变 
化 的 响应 存在 时 间 差 异 。 但 荒漠 植被 NPP 与 气温 


的 滑动 相关 系数 保持 为 负 相 关 , 与 降水 .干燥 度 的 滑 
动 相关 系数 保持 为 正 相 关 , 与 太阳 辐射 的 相关 性 不 
明显 。 也 就 是 说 ,降水 量 的 多 少 是 影响 WPP 变化 最 
主要 的 气候 因子 ,气温 对 NPP 变化 的 影响 较 小 5 。 
在 干旱 地 区 ,气温 上 升 不 仅 影响 植被 的 光合 作用 和 
呼吸 作用 ,还 加 强 了 莹 散 并 会 降低 土壤 湿度 ,限制 植 
物 的 生长 ,从 而 会 导致 NPP HFE, SAREA 
相 比 ,降水 量 的 增加 会 使 土壤 水 分 增加 ,更 有 利于 匾 
漠 植 被 的 生长 。 而 总 体 上 ,荒漠 植被 与 水 热 因 子 的 
相关 关系 在 研究 时 段 均 有 进一步 减弱 的 态势 , 即 荒 
漠 植 被 NPP 对 气候 因子 的 变化 愈 来 愈 不 敏感 。 

本 文 只 探讨 了 4 种 气候 因子 与 荒漠 植被 的 相关 
关系 ,有 一 定 的 局 限 性 。 今 后 需 加 强 车 漠 植 被 NPP 
与 不 同 深度 土壤 湿度 .土壤 特征 等 因素 之 间 关 系 的 
分 析 , 因 为 土壤 水 分 是 干旱 区 植被 生存 及 生长 发 育 
的 关键 限制 因子 。 另 外 ,人 类 活动 会 导致 荒漠 植 
被 覆盖 地 区 的 土地 利用 方式 发 生变 化 ,也 对 植被 
NPP 产生 影响 。 本 文 西北 地 区 荒漠 植被 分 布 区 是 
从 1:4 000 000 的 《中 国 植被 图 》 中 提取 出 的 ,相对 
于 34 a 的 研究 时 段 来 说 过 于 粗糙 ,今后 应 采用 多 时 
相 的 遥感 图 像 提 取 植 被 分 布 图 ,提高 研究 的 准确 度 。 
因此 ,结合 NPP 的 遥感 估算 ,充分 利用 多 时 相 \ 不 同 
来 源 的 遥感 资料 进行 植被 的 动态 监测 是 今后 要 做 的 
工作 。 
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Dynamic State of Desert Vegetation Productivity and Its Relationship 
with Water-Heat Factors in China 


LIU Xue-jia', DONG Lu’, ZHAO Jie’, DU Zi-qiang’, ZHANG Hong' 
(1. College of Environment and Resources , Shanxi University , Taiyuan 030006 , Shanxi , China ; 
2. Institute of Loess Plateau Research ,Shanxi University , Taiyuan 030006 , Shanxi , China ; 
3. College of Natural Resources and Environment , Northwest A&F University , Yangling 712100 , Shaanxi , China ) 


Abstract: In this paper, the Carnegie-Ames-Stanford Approach (CASA) model was used to estimate the NPP of 
desert vegetation in its growing season from 1982 to 2015 ,and the linear regression and GIS spatial analysis methods 
were applied to analyze the spatiotemporal variation of NPP. The purposes of the study were to explore the net pri- 
mary productivity of desert vegetation and its correlation with water-heat factors in China’s arid regions in recent 
three decades. The sliding correlation coefficient was used to analyze the relationship between NPP and water-heat 


-? «a7! and the overall level of 


factors. The results showed that: Q) The average NPP per unit area was 42 g + m 
NPP was low. Spatially ,it was higher in the northwestern and eastern margins but lower in the central , southern and 
mid-eastern regions; @ The average annual NPP of desert vegetation was 5.78 x10" g + a~'. For the interannual 
variation ,the linear change rate of the total NPP of desert vegetation was 1.64 x10" g + (10a) ~' (P =0.054) in 
China from 1982 to 2015 , desert vegetation was generally improved , but the total change trend was staged. The total 
amount of NPP showed an extremely significant growth trend (1.25 x10” g + a7',P <0.01) during the period of 
1982 — 1993 ,it was significantly decreased from 1993 to 2006 ( -6.42 x10" g - a/',P <0.01) , and increased 
slowly during the period of 2006 -2015 (1.70 x10" g +a ',P>0.05). From the perspective of spatial change , 
47.65% NPP of desert vegetation was in an increase trend, which was mainly distributed in the Alashan Plateau, 
north piedmont of the Tianshan Mountains, western marginal zone of the Tarim Basin, southeastern marginal zone of 
the Qaidam Basin,south piedmont of the Altun Mountains and Kunlun Mountains; @) From the perspective of the 
change of correlation between desert vegetation NPP and climatic factors with the time , the sliding correlation coeffi- 
cient between NPP and temperature was negatively correlated with the time, but positively correlated between NPP 
and precipitation and dryness. There was no significant variation of sliding correlation coefficient between NPP and 
total solar radiation with the time. In general, the correlation between desert vegetation and water-heat factors was 
further weakened during the study period, that is, the desert vegetation NPP was increasingly insensitive to the 
change of climatic factors. 


Key words; net primary productivity; desert vegetation; climatic factor; sliding correlation; hydrothermic factor 


